Switching between nematic filled states: flexoelectric and order-electric couplings 



M. L. Blow and M. M. Telo da Gama 

Centro de Fisica Teorica e Computacional, 
Avenida Professor Gama Pinto 2, 
P-1649-003 Lisboa, Portugal, 
and 

Departamento de Fisica, 
Faculdade de Ciencias da Universidade de Lisboa, 
P-1749-016 Lisboa, Portugal 
(Dated: August 29, 2012) 

Using numerical simulations, we model a nematic liquid crystal, in coexistence with its isotropic phase, filling 
a substrate patterned with rectangular grooves. We analyse its response to an externally-applied electric field, 
under flexoelectric and order-electric coupling, with an emphasis on identifying switching transitions between 
bistable morphological states. We identify that order-electricity provides a means of switching between states 
where the plateaux between grooves are dry and where they are wet by a nematic layer, without affecting the 
configuration of the nematic within the groove. We identify that flexoelectricity is a viable means of changing 
the topological configuration of the nematic in the groove, provided that the coupling differentiates between the 
type of distortion at the corners of the substrate. We determine precisely the field magnitudes and orientations 
required to effect each type of transition. 



I. INTRODUCTION 



The light-modulating properties of liquid crystals, arising from their anisotropic electric permittivity, together with the 
capability to change those properties by altering the orientational structure of the liquid crystal, for example by the application 
of external electric fields, make them of great importance for display devices. As well as the conventional monostable devices, 
where one display state exists in the absence of any electric field, and the others persist only in its presence, there is great 
interest in bistable devices, where two distinct display states exist in the absence of external fields, which are required only to 
switch between them, allowing greatly enhanced energy efficiency and portability. Such devices commonly exploit the multiple 
topologies afforded by structured substrates. 

The anisotropic part of the permittivity arises naturally from the tensorial nature of the apolar nematic order. This same 
apolar property precludes a nematic from possessing a bulk ferroelectric moment. However, spatial variations in this apolar 
order can themselves be polar, and thus an associated spontaneous polarisation is permitted. For example, distortions in the 
arrangement of molecules may cause them to align in accordance with a shape asymmetry. If such an asymmetry of the 
molecules is accompanied by a permanent dipole moment, the result is a net polarisation in regions of distortion. This effect 
was first identified by Meyer [T| and is termed flexoelectricity. Similarly, gradients in the degree of ordering may induce a 
polarisation, for example at defects and interfaces. This was first noted by Prost and Marcerou[2 1, and is termed order-electricity. 

There have been numerous theoretical and simulational studies of field-induced switching. Many [3-7 1 have focused on the 
zenithal bistable nematic device (ZBND), which utilises a sinusoidal substrate to achieve bistability |8], whilst others have 
considered flat substrates where the anchoring is spatially varving||9l [Toll or inherently bistable lfTTI . In all of these studies, the 
system temperature is below coexistence, deep in the nematic phase. 

Structured substrates are also of current research interest due to their wetting properties. In particular, the control and 
manipulation of interfaces is utilised in a wide range of microfluidic devices and smart surfaces. In addition to being completely 
dry or completely wet, a structure may be "filled", with an interface pinned on the topography |[T2l[T3l . Controllable undulation 
of this interface can be exploited for optical applications |[l4l and structural assembly |[T5l . When the wetting fluid is nematic, 
the behaviour is additionally enriched by the long-range interactions between the substrate and the interface, which are 
mediated by the elasticity of the nematic's internal structure, and by orientation-dependent effects (anchoring) at the interface it- 
self fl6jfT8 1|. On such surfaces, a bistability of filled states is achievable, with interfaces characterised by differing shape profiles. 



In this paper, we examine how we may switch between filled states on a structured surface with the use of electric fields. We 
model the system using Landau-de-Gennes theory with Beris-Edwards dynamics, numerically solving the equations of motion 
with the hybrid lattice-Boltzmann/finite-difference method. The paper is arranged as follows. In Sec. [II] we describe the statics 
and dynamics of a model nematic in an externally-applied electric field and in Sec. Ill we charaterise the system and the 
simulation parameters. In Sec. IV we present qualitative and quantitative results for the switching transitions between filled 
states driven by order-electric and flexoelectric couplings. We first show, in Sec. IV A that order-electricity may be used to wet 
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and dewet parts of the substrate. Such switching is equivalent to electrowetting, which is extensively documented for isotropic 
fluids lH9ll , Then, in Sec. |IV B| we consider how fiexoelectricity may facilitate switching between filled states with different 
topological arrangements of the orientational order. We provide a quantitative analysis of the transitions, identifying the range 
of applied fields that will produce them, and the free-energy barriers and timescales involved. We conclude in Sec. [V] with a 
discussion and perspectives for future work. 

II. THE MODEL 

In order to model phase separation, interfaces and defects, the required nematic order parameter is a traceless, symmetric 
tensor EOl . 

Qa/3 = \S (3n a np - 5 a p) + \B (mam/3 - l a lp) , (1) 
where {n, m, 1} is an orthonormal set, S is the degree of nematic order along n, and B is the biaxial order along m. 

The Landau-de-Gennes free energy functional, defined over the fluid domain T> and its boundary with the substrate &D is 
given by, 



J~ — / (/bulk + /elast + /diel + /flexo) d 3 r + / / wet t2 r. 
Jv JdV 



(2) 



where 



/bulk — 1^4 |t*Q Q( 3Q^ q — |(2 + T*)Q a pQf) 1 Q ia + | [Q a pQpa] 2 ^ , (3) 

/elast = ^LxdjQafjdjQafj + \L 2 d a Q al dpQ 1 p + \LzdpQa~fdaQ~fP, (4) 
/diet = — 2 e ° ( eT< W + txQap) E a Ep, (5) 
/flexo = —^XoEadpQap — |x+-Eq<9/3 (Q al Q ' 7 p) — \\'lE 1 d 1 (Q a pQp a ) ~ gX-^a (QaydfiQ^p — Q 7 pdpQ aJ ) , (6) 

/wet = —\wxQ a pQ a p + \w2QapQp0f (7) 

A is a positive constant and T* is a parameter related to temperature, which is set to 1, corresponding to coexistence between 
an isotropic and a nematic phase with uniaxial order parameters 5; so = and S nem — 1 respectively. 

The Li are related to the Frank elastic constants by 

K 1 = K 3 = l(2L 1 +L 2 + L 3 ), K 2 = \Lx. (8) 

In addition to imposing a free energy penalty to gradients in orientation, the elastic terms also penalise gradients in order 
parameter. At interfaces between the nematic and isotropic phases, the balance of elastic and bulk contributions leads to a 
characteristic interfacial width (correlation length) of 



3£i + 2(L 2 + L 3 ) 

A ' K 1 

and surface tensions, 

7 _L = i^dLx+^+La), 7|| = ^y/M^i + L 2 + L 3 ) (10) 

for the cases of homeotropic and planar alignment respectively. Comparing the surface tensions, it is seen that the preference 
of alignment, or anchoring, is dependent on the sign of L 2 + L 3 relative to L\. 

/diei is the dielectric contribution to the free energy, consisting of the usual isotropic part, with relative permittivity ej, and a 
deviatoric part proportional to Q. /fl exo relates to fiexoelectricity and order-electricity, and comprises four independent terms, 
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when taking gradients of Q up to quadratic order. It is useful to compare the polarisation described in /fl exo to a phenomenological 
description in the uniaxial framework ETI . 

P a = e + Sdp {n a np) + e_S {riadprip - npdpn a ) + e Q d a S + e A (2n a np - 5 a p)dpS (11) 

e± are the flexoelectric coefficients, corresponding to polarisations arising from the respective alignment and opposition of 
splay and bend [22|. eo and e^ are the order-electric coeffients l2TI . The eo term purely couples to the gradient of S, describing 
a polarisation across the interface independent of the alignment. The e^ term is alignment dependent, such that homoetropic 
and planar configurations have opposite polarities. 

Substition of (Eqn. [TJ into (Eqn. [6j, and comparison with (Eqn. Hi, yields GTl 



e+ = 2 Xo + Sx+ 
e_ = Sx- 



(12) 



e = 3XO + 3SX+ + SX2 

ea = Xo + Sx+ 

A set of {x} can be found to satisfy a general choice of {e} for most fixed S. But {e} cannot be consistently determined when 



S is varying; an important point of consideration given that e and ca act only on gradients in S. Eqns. ( 12 1 indicate that X- 
and xi mav be used to independently vary e_ and eo respectively, but it is impossible to set ca without affecting at least one 
of eo and e+. For example, to set a non-zero eA, while maintaining eo = across the interface, it is necessary to set xo = 
and 5%+ = —3x2, in which case e+ does not vanish in the bulk. Alternatively, e + may be eliminated in the bulk by setting 
X+ = — 2xq> but this does not allow for eo to be zero everywhere across the interface. To remove these interdependencies, it is 
necessary to include higher order terms in the polarisation, but this shall not be considered here. The restriction may be likened 
to the fixed equality between K\ and if 3 expressed in Eqn. [8] 

Finally, w\ and w 2 are the anchoring strengths, and Q s is the preferred form of the Q-tensor at the substrate. 

Minimisation of the free energy in Eqn.Q, subject to the constraints that Q be symmetric and traceless, leads to the equilib- 
rium condition H vo1 = throughout V, and H sulf = along dT>, where 



H-ofi = |^ {GQa-yQ-yP ~ 2Q-ysQs~/daf3 ~ Qafi ~ ^Q-ySQs-yQaf}) + L\d 11 Q a p + (L 2 + L^d^gP af3~/S 

+ \^of-T>{.E a Ep — ^8 a pE 2 ) — |xo (^d a E/3 + \dpEa) — ^S a pd 7 E 7 ) 

- 1x2Qaj3d 1 E 1 - |x+P Q /3 7 5 {d 1 E s + d s E 7 ) - |x_ {E s d 1 - E 7 d s ) P a/ 3 7 s (13) 

and, 

H af3 = -2Li^ 7 a 7 Q Q/3 - (L^dg + L 3 V S d 7 ) P Q/ 3 7 <5 + f Wi<2^3 + I {X2E 1 V 1 - W 2 ) Qaf3 

+ f Xo {\E a vp + \vpE a - l5 al3 E~ ( v~ ( ) + P Q( 3 7 5 (|x+ {v~/E s + E 7 v s } + |x- {^E s - E~ ( v s }) (14) 
where P Q ^ 7 a = \8 ai Qsp + \Q a s& 1 (3 — \5 a pQ 7 s and v is the inward normal to the substrate. 

Of note is that, in the bulk, only the sum L 2 + £3 is of consequence. The distinction between L 2 and L3 only manifests itself 
in the surface terms, where it fine-tunes the line tension. Similarly, in the bulk, the x- coupling acts only via gradients in Q 
(using curlE = 0), while \+, Xo, Xi act only through gradients in E J7). In a uniform field, these latter couplings produce only 
surface effects. 

To describe the dynamics of the fluid, we must introduce density and velocity fields p and u. The evolution of the fluid with 
time t is described by the continuty, Navier-Stokes, and Beris-Edwards equations. 

d t p + dp(pu ) = Q (15) 
p (d t + Ufi d ) u a = d p [p{2 v A aP - 5 a pc 2 s } + {{Z afhS + T afhS } H;f] + H%d a Q 7fl (16) 
(d t + u 7 <9 7 ) Q a p = —(EapysAjs — T Q( g 7 (5f2 7 5 + TEP°p (17) 

where 

T Q) g 7 a = Q ai 5ps — Sa-yQfss, = 2 {dpu a + d a up) , f2 Q(3 = \ (dpu a — d a up) . 



(18) 
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FIG. 1: (a) States where the nematic phase fills the crenel: (i) BB, (ii) BS, and (iii) SS. (b) States where the nematic phase fills the crenel and 
wets the plateau, but is pinned at the upper comers of the groove: (i) BB, (ii) BS, (iii) SS. 



r\ is the isotropic part of the kinematic viscosity (there are also anisotropic contributions to the viscosity from the interplay 
between Eqns. ( To*p7 i, c s is the speed of sound, £ is a parameter depending on the molecular details of the liquid crystal and T 
is the mobility. The system is closed by the boundary conditions (5 a p — v a vp)up = (non-slip) and i?™ f = 0. 



III. THE SYSTEM 



We consider a system confined between two neutrally wetting walls (wi = ui2 — 0) that lie perpendicular to the y direction 
with separation <C 1. As such, the system is considered to be two dimensional, lying in the xz plane. In the x direction, the 
system is periodic, whilst in the z direction, it is confined between two substrates. The upper substrate is flat and neutrally 
wetting, but the lower substrate is patterned and has a finite anchoring strength. Following [18|, we consider a substrate that 
is crenellated with rectangular grooves. Each of pair of grooves are separated by a plateau of width 20£, and the grooves have 
width 20£. The height/depth of the grooves is 10£, and the vertical distance from the plateaux to the upper substrate is a further 

We use the simulation parameters A = and L\ = Jj, which through Eqn. J9b produce a £ of 3 lattice spacings. Again 
following [18 1, we choose L 2 = 2Li and L 3 — 0. This choice of elastic coefficients gives rise to planar anchoring at the 
interface. At the substrate, we choose = § {y a vp — 5^/3), w\ = j^£ _1 ^i EH, and u>2 = 0, so that the substrate favours 
the nematic phase over the isotropic phase, and homeotropic alignment over planar. 

We focus our attention on states where the base and sides of the grooves are wetted by the nematic phase, but the plateaux 
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and upper surface are in contact with the isotropic fluid. Thus, the interface connects the two upper corners of the groove. We 
refer to these as "filled" states. Within the groove, the nematic will be distorted due to conflict between the anchoring of the 
three walls and the interface. At each of the bottom corners, the homeotropic anchoring to the rightangled pair of walls causes 
the director to either bend or splay. Applying these possibilities to the two lower corners, there are four possible filled states, 
which we label BB, SB, BS and SS, where the first and second letters indicate the type of distortion in the left and right lower 
corners respectively. Of course, SB and BS are degenerate, being mirror images of one another. 

The distortion-mediated interactions between the substrate and nematic-isotropic interface determine the shape of that 
interface and the relative free energies of the filled states. In the BB state, the interface follows the contours of the bend, 
as illustrated in Fig. [T^(i), whilst in SB/BS (Fig. [lj(ii)) the interface is flat. In the terminology of [18|, SS and SB/BS are 
thus referred to as the bent and unbent states (referring to the shape of the interface), but BB is not mentioned, because it is 
not a global free energy minimum for any of the parameters that the authors consider. Fig. [T^(iii) offers an explanation: the 
configuration in the nematic forces the alignment at the interface to be homeotropic, thereby increasing the surface tension. 

For the parameters specified for this system, lfl8l show that SS is the state with lowest free energy. We further find that 
SB/BS has a lower free energy than BB. 

There are other states in addition to the filled states. For example it is possible for the fluid region of the system to be 
completely isotropic (dry state) or completely nematic. In ||T8l . where an open system is considered, the state possesses an 
interface far from the substrate. In the system considered here, the nematic wets the upper substrate, and is thus confined in a 
finite space, so the state is not "wet" in the usually understood definition. As with the filled states, there are multiple types of 
such a state, depending on the nature of the distortion at the four corners. However, since this paper focuses on filled states, we 
do not distinguish the "wet" states in detail, and simply refer to them collectively as W. 

Another possible type of state, as identified by ifTSl . involves the entire lower substrate being wetted by nematic, but the 
interface remaining pinned at the upper corners of the groove. Such a state possesses the same morphology as the filled states 
plus a thin film on the plateau, as illustrated in Figs.[TJ). Hence we denote such states by BB, BS, SB, SS. 



To simulate the system, we use a hybrid method, solving Eqns. ( jT3] > and (16 1 with a Lattice Boltzmann algorithm, and 



Eqn. ( fTT) by a finite difference algorithm. These algorithms are evaluatated in tandem, with an additional predictor-corrector 
step to improve stability, as detailed in [25]. Most Lattice Boltzmann studies have been carried out with T* deep in the nematic 
phase[5, 9|[24][25). Those dealing with coexisting nematic and isotropic fluids have looked either two-substance systems, with 
the relative amounts of each type of material fixed E6l . or single-substance systems at or close to T* = 1 [ 27 1 , so that nematic 
and isotropic coexist in variable proportion. The system we consider is of the latter type. 

Viable simulation parameters must be mapped to physical quantities by an appropriate choice of scaling. Suitable real values 
of parameters are on the order, A ~ 10 4 Nm" 2 and 3Li + 2L 2 ~ KT 10 N, leading to f - 10~ 7 m. By comparing this physical 
value to the value calculated in simulation units above, we can determine the lengthscale corresponding to a lattice spacing as 
being ~ 10 _ 'm. 

To ensure stability, the simulation density p is limited to a minimum value ~ 1. Using the physical value p ~ 10 3 kgm -3 , 
we map a simulation mass unit to 10 _18 kg. The timescale may be determined by mapping either surface tension or viscosity. 



Using Eqn. ( 10 1, we derive a value 7 = 0.022 in the simulation units and 7 ~ 10 3 Nm 1 in physical units. These are matched 



by one iteration corresponding to ~ 10 s. 

Using a relaxation time r = 2 in the Lattice Boltzmann algorithm provides an isotropic kinematic viscosity of 0.5, which 
corresponds to 77 ~ 10 _6 m 2 s _1 . This is an order of magnitude smaller than accepted values (generally, stability restrictions 
make it difficult to provide a compatible viscosity and surface tension, without resorting to a prohibitively large grid |26|). We 
set T = 0.45 and £ = 0.4, providing deviatoric components of viscosity of a like magnitude. 

Finally, we consider the electrical parameters of the system. Typically €\ — 10, and for simplicity we assume that = 0. 
The flexoelectric coefficients are typically \ ~ 10 _11 Cm _1 [28, 29|. Overall, the dielectric and flexo/order-electric effects 
depend on eE 2 and \E respectively, so changing E changes their relative strengths. We assume that the isotropic dielectricity 
dominates, thus the flexoelectric polarisation does not affect the applied E-field, and it is not required to solve Laplace's equation 
at every step of the simulation. 




IV. SWITCHING TRANSITIONS 



A. Polarisation at the interface 



We first consider a system where only \2 is non-zero, corresponding to a system where there is a polarisation across gradients 
of S, independent of director alignment. The resulting polarisations for the BB (stable), BS (metastable), and SS (metastable), 
filled states are shown in the third column of Fig. [2] 

Figs.[3^-c show the effect of applying a field X2E = (0, —0.03). Starting in BB, a homeotropic wetting layer grows on the 
upper surface of the plateau, until the system resembles BB. The morphology within the crenel is little changed. When the field 
is switched off, the system remains in BB, this being a metastable state at the chosen parameters. Thus, we have effected a 
transition BB ->• BB. 



The behaviour can be explained as follows. As indicated in Eqn. ( 14 1, the \2 gives rise to a surface term with the same form 
as the quadratic wetting potential, with an effective W2 given by — X2"-E, and thus dependent on the orientation of the substrate 
relative to the field. On the horizontal sections of the substrate, there is a contribution in QapQpa, which is negative when the 
substrate faces upwards and hence favours S = 1, thereby causing the growth of a nematic layer on the plateau. On the vertical 
sections, however, there is no change to the effective wetting potential, so the contact line does not move. Furthermore, because 
the field is uniform, there is no force on the polarised interface. 

By reversing the direction of the field, the sign of the effective W2 changes, favouring the isotropic phase S = at the plateau. 
Figs. [3ji-f shows how applying X2E = (0, 0.03) dewets the plateau. When the field is switched off, the cycle is completed, with 
the system moving back to the starting configuration of Fig. [3^. 



We observe that transitions between BS and BS, and between SS and SS may be brought about in in a similar way (not shown). 
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We now consider how we might switch the system between the topologies BB, BS and SS. For example, the transition 
BB — > BS requires symmetry breaking, indicating that a horizontal field might be used. 

Fig.|4]shows the results of applying X2E = (0.03, 0). By the same token as above, the lefthand wall of the crenel now favours 
S = 0, while the righthand wall has (additional) preference for 5 = 1. The effect is thus that the lefthand wall dewets, causing 
the interface to sit diagonally in the crenel. However, when the field is switched off, the nematic rewets the wall, returning to 
BB without any transition taking place. 

Examining the director orienation in Fig.|4]illustrates why this approach will not cause a switch to BS. In each lower corner, 
the director remains in a bent configuration. It is unsurprising that coupling via \2 will not bring about a change in the director 
topology, since it interacts only with the degree of ordering S. 



B. Polarisation in the bulk nematic 



We have seen that pure order-electric coupling through the coefficient \2 — isn't capable of changing the director topology 
at the comers of the crenel, and hence bringing about transitions between BB, BS and SS. We thus turn our attention to the 
other coefficients present in Eqn. d6). A comparision is shown in Fig. [2] 

Examining the polarisation in the lower corners of the crenel, we note that for xo an d X+ coupling, the polarisation always 
points outwards at the corners, regardless of whether the director morphology is splayed or bent, with the result that the 
polarisation fields of the three states bear little difference from one another. We therefore conclude that coupling of this nature 
cannot strongly drive transtions between the states. 

However, when coupling is of the \- form, we see that the polarisation points outwards at corners where the director is bent, 
but inwards where the director is splayed. As such, the three states possess different polarisation fields, such that the application 
of an electric field will favour one over the others. We therefore take \- to be the relevent coefficient in driving transitions 
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FIG. 5: (a-h) The transition BB — > BS, through the application of the field X-E = (0.0192, —0.0111) and (h-l,a) The reverse transition 
BS -» BB, through X-E = (-0.0192, 0.0111) 



between the states, and set the other \ t0 zero. 



As indicated in Eqn. 13 this is the only flexoelectric coefficient that acts in the bulk when the electric field is uniform. The re- 
sulting polarisations for the bent BB (stable) and unbent BS (metastable) filled states are shown in Fig.|2j5. In BB, shown in c(i), 
the polarisation points outwards from both lower corners, whilst in BS, the polarisation in the bottom-left corner is also pointing 
outwards, but that in the bottom-right corner is instead pointing inwards. Thus the application of E in the downward-right di- 
rection would be expected to thermodynamically favour BS over BB, and may offer a means of effecting the transition BB — > BS. 

Figs.|5^-f show the effect of applying a field X-E = (0.0192, —0.0111), starting with BB in Fig. [5^. The polarisation in 
the bottom-right corner opposes the field, so this region of the nematic carries a free energy cost. The interface is thus pulled 
towards the corner in order to shrink this region, as shown in Fig. |5j). This happens until the righthand wall of the crenel has 
been dewetted, as shown in Fig. [5}:. At the bottom-right corner, the director rotates anticlockwise, passing the vertical position. 
The nematic then rewets the wall, as shown in Figs.[5}i,e, but this time creates a region in which the polarisation is aligned with 
the field. The system comes to rest in the configuration shown in Fig. [3]F, which is topologically the same as BS, but with a 
bulging interface. When the field is swiched off, the system relaxes to BS, as shown in Fig. |5^, thus completing the transition 
BB^BS. 



In order for the transition to occur, the field strength must be within the correct range. If X-E 1S t°° small, then the distortion 
of the interface will not be great enough to allow the director to rotate. For example, switching on the field may produce a state 
similar to Fig.[5|), which will then relax back to BB when the field is removed. Alternatively, if X-E is to ° large, then the entire 
system will wet. This usually occurs by the bulge of the interface in Fig. [3J being sufficient for the interface to depin at the 



upper corners. As discussed in Sec. Ill the distorted nemtic has a lower free energy density than the isotropic phase, so wetting 



of the entire system readily happens once the interface depins. 

In addition to being a suitable strength, the field must be directed in an appropriate angle. In Fig. [5] an angle of 30° below 
the horizontal is used, which offers a large range of applicable field strengths, but the transition is possible from 22° above the 
horizontal to 68° below. Fig.|6^(i) shows the outcome of applying a field with components x~(E x , E z ) as marked on the axes. 
In each case, the field is applied, and the system is left to relax to its new equilibrium. Then the field is removed, and the system 
is once again left to relax, with final state being recorded. The transition BB — > SB of course occurs under the same criteria as 
BB — > BS, but with the x component of the field reversed. 

Fig. |6^(ii) shows plots of the free energy (excluding the electrical contributions) of the system against time, for the 
fields corresponding to the respective coloured markers on Fig. |6^(i). The E field is kept on until t = 80, 000 iterations, 
within which the system is fully relaxed in all cases, although in many cases, only a much shorter time is required for the 
transition. After this time, there is no field. The red curve corresponds to parameters where the transition does not proceed, 
and the free energy returns to its original value. The brown curve indicates a situation where the system wets. We note 
that the wet state actually has the lowest free energy, due to the interface being removed, but it also corresponds to the 
highest free energy barrier. All other curves correspond to the transition BB — s- BS, hence the final free energy is higher 
than the original, because BS is metastable for this system. The blue curve shows the evolution of the system presented in Fig.{5]l. 



Having outlined the transition BB — > BS, we now consider the reverse, BS — > BB. As the down-right field favoured 
the polarisation of BS over BB, an up-left field should favour the reverse. Figs. |5^-k show the effect of applying a field 
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FIG. 6: (i) The resulting state following the application of an electric field, plotted against \-{E x , E z ), starting from the state (a) BB, (b) 
BS and (c) SS. (ii) Free energy (excluding the electrical contributions) against time, for the fields corresponding to the respective coloured 
markers in (i). The field is switched on until t — 80, 000 iterations, after which time it is switched off. The blue curves (inverted triangle) in 
(a) and (b) correspond to the evolution depicted in Fig. [5] The purple curve (unfilled circle) in (b) corresponds to Fig. [7] The blue (inverted 
triangle) and green (diamond) curves in (c) represent Figs.[8]and|9]respectively. 
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FIG. 7: The transition BS — > SS, through the application of the field X-E = (0, —0.0315). The field is switched on in snapshots (a-d) and 
off in (e-f) 
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FIG. 8: The transition SS — > BB, through the application of the field X-E 
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FIG. 9: The transition SS -> BS, through the application of the field X-E = (0.000193, 0.00221). The field is switched on in snapshots (a-d) 
and off in (e-f) 

E = (—0.0192, 0.0111) to the state BS. Its effect is not to pull the interface, as was the case in the previous transition, but to 
rotate the director until the geometry is such that the polarisation is aligned to the field. Once the director at the right-hand 
wall has got past being vertical (as seen in Fig.. the system readily evolves until it is in a state topologically similar to BB, 
as shown in Fig. [3], but with the interface higher and relatively unbent. Switching off the field causes the interface descend 
(Fig.|5jc,l,a) until it has returned to the bent configuration. 

Fig. [6j3(i) shows the field response of the system, when initially in state BS. We note that the transition BS — > BB is effected 
by applying the electric field in the upward and leftward directions, unsurprisingly the opposite of BB — > BS, and with a range 
of fields similar in order of magnitude, if slightly smaller. It is not possible to switch from BS to its mirror image SB directly 
through a single field application. However, it can be done through two successive field applications and relaxations, with BB 
as an intermediate state. 

We notice also a small region in the bottom-left corner of Fig. |6j)(i), in which the transition BS — » SS occurs. Fig.|7]outlines 
the mechanism. Just as an upward field caused the director to rotate on the righthand wall, Figs. [5^-c show that a downward 
field causes a rotation of the director on the lefthand wall, turning the bent topology in the bottom-left corner into splay. The 
interface then bulges upward, as shown in Fig. [5}l, and we note that wetting will occur if it reaches the ceiling of the system. 
Indeed, if the ceiling were higher, then the BS — > SS region of Fig. ^>(i) would extend further. When the field is swiched off, 
this bulge relaxes, as shown in Figs.|5^-f. 

Fig. |6j3(ii) shows plots of the free energy (excluding the electrical contributions) of the system against time, for the fields 
corresponding to the respective coloured markers in Fig. |6j3(i). As with Fig. |6^(i), the E field is kept on until t = 80, 000 
iterations, and is then switched off. The red curve corresponds to a case where no transition occured, the orange, green and blue 
curves to various instances of BS — )• BB, the purple curve to BS — > SS, and the brown curve to a scenario where the system wets. 

Finally, we investigate the response of the state SS to a field, which is shown in Fig. [6j;(i). Immediately of note is that, if 
the field has an upward component, only a very small magnitude is required to enact a transition, either to BB if the field is 
close to vertical, or to BS or SB if it has a sufficient lateral component. Only against downwardly-orientated fields is SS robust. 
Fig.|6};(ii) shows the free energy plots. As indicated by the green and blue curves, which correspond to points very close to the 
origin, there is only a very small energy barrier between SS and either BB or BS/SB. 

V. DISCUSSION 

We have used a hybrid lattice-Boltzmann/finite-element approach to simulate a nematic fluid filling grooves in a confined 
rectangular structure, and have discovered that the application of electric fields, with flexoelectric and/or order-electric coupling 
to the nematic, may be used to switch between morphologies. We have found that order-electric coupling enables switching 
between states where the plateaux of the structure are dry and those where the plateaux are wet, with the structure in the groove 
being unchanged, while flexoelectric coupling, specifically the component with coefficient allows for switching between 
the various topologies that may fill the groove. 
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At the end of Sec. Ill 



we noted that we apply the assumption of uniform E field throughout the system. This greatly 
simplified the simulation procedure, by allowing field gradient terms in Eqn. [13] to be neglected, and removing the need for 
the solution of electrostatic equations. We return to consider the viability and limits of a such an assumption, by considering 
scenarios where it might not hold. 



Firstly, if the substrate permittivity is different to ej, then a non-flat substrate will distort the electric field around it. Such an 
effect is considered in many other works on flexoelctric switching j5]-[7l, and in some cases plays an essential role in switching. 
However, it is usually possible to choose materials so that the difference in permittivities is small, and as we have shown, 
switching can be achieved without exploiting this effect. 

Another potential complication arises from the effect of the polarisation itself, which would lead to E being not only nonuni- 
form, but changing with the nematic configuration. Maxwell's equations in media give 



d a {e (eiS a fi + e A Q a p) Ep ~ ^xodpQap - \x+dp (Qa-yQ-yp) ~ |x2<9 Q {Qp 1 Q 1 p) - %X- {Qa-ydpQ-fP - Q 7 pdpQ a ~,)} = 
£o {tid a E a + e A Q a pd a Ep) = —coeAEpd a Q a p + ^Xod a pQ a p + ^X+d a p (QayQjp) + \x2daa (QptQ^p) 

+ §x_<9 Q {Q ai dpQ 1 p - Q 1 pdpQ ai ) . (19) 



Solution of Eqn. ( fT9] l at every timestep would add significant computational burden, not least because, unlike Eqns. ( fl"5fT7 1, it 



is not local and would require repeated iterations over the grid. There has only been one lattice Boltzmann study where this has 
been done [24|. However, it is consistent with an approximately uniform field provided that €\ eo and e eiE x£nem£ _1 - 
Whether the first condition is met depends solely on the material properties (for example, [20 1 quotes a ratio of eo/ej ~ 0.1 



for MBBA), while the second condition also depends on the field strength. In Sec. IV we ascertained that the transitions are 
effected when X-E ~ 10 _3 Nm~ 1 in physical units. Assuming that the first condition is met, the second may be expressed as 
the ratio, 



^ X— 'S'nem (20) 

£e eiE 

A typical value X- ~ 10 _11 Cm _1 ll28l|29l , requiring E ~ 10 8 Vm _1 and resulting in r ~ 0.01. This value o r validates the 
second condition, while the potential required to produce E on the scale of the device, is on the order of volts, which is feasible 
for many applications. 



We have studied the limit where the channel is confined in the y direction, and the system is thus two-dimensional. It is 
interesting to consider the opposite limit, where the channel extends indefinetly in the y direction. In the case of zero field, 
translational symmetry along y dictates that the states are unchanged from the narrow limit |18|. However, the application of 
an electric field may induce gradients in Q along y, breaking the translational symmetry to produce periodic variation in this 
direction. One difficulty in simulating such a scenario lies in setting the length of the simulation box along y, since this sets the 
periodicity of any variations along y, and may not match that naturally taken by the systems. Calculations may give an indication 
to this repeat length [7 ], and in terms of identifying equilibrium states, it is possible to identify it by varying the length of the sim- 
ulation box and finding minima in the free energy. However, in a dynamical situation, the natural periodicity may not be constant. 

A related three-dimensional system to investigate, which would not entail such difficulties, would be a substrate patterned 
with cuboidal holes, i.e. comparable spacings on the x and y direction. Such a system would have radically different states, both 
in the absence and presence of an electric field. 
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